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ABSTRACT 



33 We report the discovery with Fermi /LAT of 7— ray emission from three radio- 
loud narrow-line Seyfert 1 galaxies: PKS 1502+036 (z = 0.409), 1H 0323+342 
(z = 0.061) and PKS 2004-447 (z = 0.24). In addition to PMN J0948+0022 
(z = 0.585), the first source of this type to be detected in 7 rays, they may form 
an emerging new class of 7— ray active galactic nuclei (AGN). These findings can 
have strong implications on our knowledge about relativistic jets and the unified 
model of AGN. 

34 Subject headings: quasars: general - galaxies: active - galaxies: Seyfert - gamma 

35 rays: observations 



36 1. Introduction 

37 The advent of the Fermi Gamma-ray Space Telescope (hereafter, Fermi), with its excel- 

38 lent performance, is changing our perception of the sky at high-energy 7 rays. Specifically, 

39 the recent detection in the MeV-GeV energy range of the radio-loud narrow-line Seyfert 1 

40 (RL-NLS1) quasar PMN J0948+0022 strongly supports the presence of a closely aligned 

41 relativistic jet in this peculiar system (Abdo et al. 2009a, c; Foschini et al. 2009a). This 

42 is quite surprising, since NLSls are generally hosted in spiral galaxies and the presence of 

43 a fully developed relativistic jet is contrary to the well-known paradigm that this type of 

44 system is associated to ellipticals (e.g. see Marscher 2009 for a recent review). 
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45 Already in the past, several authors inferred from the multiwavelength properties of 

46 RL-NLSls some parallelism between this type of source and blazars (e.g. Komossa et al. 

47 2006, Yuan et al. 2008, Foschini et al. 2009b), suggesting similarities with different flavors 

48 of the blazar types (quasars or BL Lacs). What was missing in these previous studies is the 

49 7— ray detection, which is important to confirm the presence of a relativistic jet, to measure 

50 its power, and to study the characteristics of this type of sources by modeling their spectral 

51 energy distributions (SEDs). Indeed, blazar-like radio emission from radio-quiet AGNs has 

52 been already reported (e.g., Brunthaler et al. 2000, Brunthaler et al. 2005, Lister et al. 

53 2009; see also the review by Ho 2008), but no 7 rays have been yet detected from these 

54 sources. This suggests some evident differences between fully developed relativistic jets (i.e. 

55 emitting over the whole electromagnetic spectrum, from radio to 7 rays) and jet-like (perhaps 

56 aborted?, cf Ghisellini et al. 2004) structures, whose emission at 7 rays in the MeV-TeV 

57 energy range - if any - has not been yet reported. 

58 Understanding these differences could give important insights into jet formation. There- 

59 fore, we started a larger program aiming at the detection of 7— ray emission from other RL- 

60 NLSls. Since no complete sample of RL-NLSls is available in the literature, we have built 

61 a sample by merging all the sources of this type in the catalog of Zhou & Wang (2002) and 

62 in the lists of Komossa et al. (2006) and Yuan et al. (2008). We have adopted a threshold 

63 in the radio-loudness Rl = GHz/,/440 nm = 50, in order to avoid doubtful sources at the 

64 border (R L = 10) with radio-quietness. The optical properties are those typical of NLSls, 

65 i.e. narrow permitted lines (FWHM R(3< 2000 km s" 1 ), [OIII]/H/3 < 3 and the bump of Fell 
ee (see Pogge 2000 for a review). At radio frequencies, RL-NLSls display strong and variable 

67 radio emission, with brightness temperature well above the inverse- Compton limit (see, e.g., 

68 Yuan et al. 2008, Table 3). 

69 This selection resulted in a list of 29 sources, most of them selected from the Sloan Digital 

70 Sky Survey (SDSS). This low number should not be a surprise. According to Komossa et 

71 al. (2006), RL-NLSls are more rare than quasars: 7% of NLSls have Rl > 10 and only 

72 2.5% exceed Rl = 50, while generally 10-20% of quasars are radio-loud. Although based 

73 on a small sample of objects (128), these results have been confirmed by studies on larger 

74 samples from the SDSS (Zhou et al. 2006, Whalen et al. 2006). It is worth noting that this 

75 sample is not complete, although the fact that it is based on the SDSS ensures a minimum 

76 degree of representativeness. 

77 Here we report three new detections at 7 rays of RL-NLSls obtained with Fermi /LAT. 

78 In addition to the already known PMN J0948+0022 (reported by Abdo et al. 2009a,c; 

79 Foschini et al. 2009a), this increases the number of 7— ray detections of NLSls (3 + 1=4 
so out of 29 sources, 14% of our sample), suggesting that they form a new class of 7— ray 
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81 emitting AGNs. 

82 Throughout this work, we adopted a ACDM cosmology from the most recent WMAP 

83 results, which give the following values for the cosmological parameters: h = 0.71, Q m = 0.27, 

84 Q\ = 0.73 and with the Hubble-Lemaitre constant Hq = lOO/i km s _1 Mpc -1 (Komatsu et 

85 al. 2009). 

86 2. Data Analysis 

87 The data from the Large Area Telescope (LAT, Atwood et al. 2009) onboard Fermi 
sa were analyzed following the same scheme described in Abdo et al. (2009a), but using a more 

89 recent version of the software (Science Tools v 9.15.2), Instrument Response Function 

90 (IRF P6_V3_DIFFUSE, Rando et al. 2009) and background subtr action (0). The data analyzed 

91 span between MJD 54682 (2008 August 4) and 55048 (2009 August 5). 

92 Four LAT sources with T£0 greater than 25, which is equivalent^ to ~ 5a, can be 

93 associated with RL-NLSls by using the Figure-of-Merit (FoM) method outlined in Abdo 

94 et al. (2009b) and already used to associate other LAT sources with radio counterparts. 

9 5 One of them is the already well-known PMN J0948+0022 (Abdo et al. 2009a). The three 

96 remaining new sources were associated to the RL-NLS1 objects PKS 1502+036 (z = 0.409), 

97 1H 0323+342 (z = 0.061) and PKS 2004-447 (z = 0.24). The results are summarized in 

98 Table ffl 

99 In order to build SEDs, we retrieved all the publicly available data of these sources. In 

100 the case of PKS 1502+036, no X-ray data were found and, therefore, we asked for a 5 ks snap- 

101 shot with the Swift satellite, which was performed on 2009 July 25 (ObsID 00031445001). We 

102 retrieved from the public Swift archives, optical/UV/X-rays data also for PMN J0948+0022 

103 (Abdo et al. 2009a) and 1H 0323+342. For the latter, we analyzed all the observations 

104 performed between 2006 July 6 and 2008 November 16 (15 observations) and averaged the 

105 results (total exposure on XRT was 57 ks). In the case of PKS 2004-447, one XMM-Newton 
we observation was found (ObsID 0200360201, performed on 2004 April 11; 42 ks exposure on 
107 EPIC, see Gallo et al. 2006) and analyzed. 



1 We adopted the same version of software and calibration database that is publicly available at 
http : //f ermi . gsf c .nasa.gov/ ssc/dat a/analysis/software/. 

2 Likelihood test statistic, see Mattox et al. (1996) for a definition of this test. 

3 The significance of the detection in a is roughly \/TS. 
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Table 1: 7-ray characteristics of RL-NLSls detected by LAT with TS > 25. 



Name 


z 


a 


6 


r 95% a 




r c 


TS 




Ass.° 






[deg] 


[deg] 


[deg] 


[10~ 8 ph cm- 2 s- 1 ] 




[ 


10 46 erg s- 1 ] 


[%] 


1H 0323+342 


0.061 


51.27 


+34.19 


0.21 


5.7 + 0.3 


2.74 + 0.03 


78 


0.02 


92 


PKS 1502+036 


0.409 


226.26 


+3.37 


0.18 


7.5 ± 0.3 


2.81 + 0.03 


186 


2.1 


93 


PKS 2004-447 


0.24 


301.99 


-44.51 


0.09 


2.3 + 0.3 


2.5 + 0.4 


43 


0.2 


81 


PMN J0948+0022* 


0.585 


147.28 


+0.39 


0.12 


14.6 ± 0.9 


2.78 + 0.06 


647 


11 


99 



"Error radius at 95% confidence level. To take properly into account systematic effects in the calculation of 
the error radii, we added the absolute systematic error of 0°.0075 in quadrature to the 68%, then we convert 
to 95% by multiplying by 1.62 and, finally, we multiply by 1 .2, which is the relative systematic error (Abdo 
et al., in preparation). 

b -f— ray flux for E > 100 MeV. We quoted only statistical errors, while the systematics to be considered are 
10% at 100 MeV, 5% at 500 MeV and 20% at 10 GeV (Rando et al. 2009). 
c Photon index of the power-law model used to fit LAT data. 
^Observed 7— ray luminosity. 

e Confidence level of the association according to the FoM. 
f See Abdo et al. (2009a). 

108 PMN J0948+0022 was extensively studied in Abdo et al. (2009a,c), while the data of 

109 1H 0323+342 and PKS 2004-447 were already discussed in Foschini et al. (2009b). Partic- 

110 ularly, we found in the case of 1H 0323+342, there is evidence of spectral variability in the 
in X-ray energy band, even on timescales of a few days. The photon index is generally soft, as 

112 for a typical NLS1 (cf Leighly 1999), but sometimes it displays a hard tail (r ~ 1.3 — 1.6) 

113 during a simultaneous increase of the optical/UV flux (Foschini et al. 2009b). 

114 PKS 1502+036 was first observed in X-rays on 2009 July 25 (Swift ObsID 00031445001). 

115 No previous X-ray observations were available and even ROSAT measured just an upper 
lie limit, which was not very stringent (< 2 x 10~ 12 erg cm -2 s -1 in the 0.1 — 2.4 keV, see 
117 Yuan et al. 2008). We based our plans for the observation taking into account the similar 
us optical characteristics shared with the better known PMN J0948+0022 (Abdo et al. 2009a). 

119 However, the effective exposure of 4.6 ks on XRT was not long enough to collect sufficient 

120 photons to apply the x 2 test statistic. Therefore, we grouped the photons to have at least 10 

121 counts per bin and applied the Cash statistic, which allows the parameters to be estimated 

122 through the likelihood ratio (Cash 1979). 

123 In order to make a homogeneous data set, we have reanalyzed all the Swift and XMM- 

124 Newton data adopting the same procedures outlined in the previous works, but using the 

125 most recent software packages version (HEAS0FT v . 6.6.3 for Swift and SAS v . 9.0.0 for 

126 XMM-Newton) and calibration databases (updated on 2009 June 5 for Swift and 2009 June 
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Table 2: Summary of the Swift and XMM-Newton analyses. 











X-rays 










Source 


Inst a 


Exp b 


iV H c 


r d 




r 2 f 


FS 


Stat' 1 


1H 0323+342 


X 


57.1 


11.7 


2.03 + 0.02 1 






14.0 + 0.1 


1.06/313 


PKS 1502+036 


X 


4.6 


3.89 


1 0+ ' 9 
1 - u -0.8 






0.43 + 0.10 


2 


PKS 2004-447 


E 


15,22,22 


2.96 


2 0+°' 2 


0.67 + 0.08 


1.50 + 0.03 


1.5 + 0.1 


1.05/315 










Optical/UV 








Source 


Inst a 


Ay 1 


V 


b 


u 


uvwl 


uvm2 


uvw2 


1H 0323+342 


U 


0.62 


15.82 + 0.03 


16.38 + 0.03 


15.50 + 0.03 


15.82 + 0.04 


16.01 + 0.04 


15.92 ± 0.04 


PKS 1502+036 


u 


0.21 


19.0 + 0.2 


19.5 + 0.1 


18.6 + 0.1 


18.54 + 0.08 


18.38 + 0.08 


18.40 + 0.06 


PKS 2004-447 


o 


0.16 




18.91 + 0.08 


18.54 + 0.09 


18.65 + 0.08 


19.0 + 0.2 





instrument used for the observation: X for XRT and U for UVOT onboard Swift; E for EPIC PN, MOS1, 
and MOS2 and O for OM onboard XMM-Newton. 

h Net exposure in kiloseconds. In the case of EPIC, the exposures of PN, MOS1 and MOS2, respectively, are 
reported. 

c Galactic absorption in units of 10 20 cm~ 2 from Kablerla et al. (2005). 

d Photon index of the power-law model or low-energy photon index in the case of broken power-law model. 

e Break energy in the case of broken power-law model [keV]. 

•^High-energy photon index in the case of broken power-law model. 

^Observed flux in the 0.2 — 10 keV energy band in units of 10~ 12 erg cm~ 2 s _1 . 

''Statistical parameters: x 2 /dof for \ 2 or number of PHA bins for the Cash statistic. 

4 It is worth noting that the addition of an unresolved gaussian emission line at E = 6.5 ± 0.3 keV with 
equivalent width <~ 147 eV gives a A% 2 = 9.3. 

'Magnitudes absorbed by the Galactic column in the visual filter (from Kalberla et al. 2005) and used as 
reference to calculate the absorption in other filters by means of the Cardelli et al. (1989) extinction law. 
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127 11 for XMM-Newton). X-ray data were fitted to power-law or broken power-law models with 

i2B Galactic absorption from Kalberla et al. (2005), while optical/UV data were dereddened for 

129 the proper values calculated according to the extinction laws by Cardelli et al. (1989). None 

130 of the sources showed evidence of X-ray absorption in excess to the Galactic value. The X-ray 

131 photon indexes are generally in agreement with those typical of blazars (c.f., for example, 

132 Foschini et al. 2006; Maraschi et al. 2008). The results are summarized in Table [2j All the 

133 results of the present reanalysis are consistent with the previous ones already reported in 

134 Foschini et al. (2009b). 



Table 3: Parameters used to model the SEDs. 



Name 


p a 


logAf b 




p/d 


L d c 


B l 


r b g 


0v h 


Tc, break 


<v 1 

/c.max 


si m 


S2 n 


1H 0323+342 


1.9 (650) 


7.0 


116 


1.0 


1.4 (0.9) 


30 


12 


3 


60 


6000 


-1 


3.1 


PKS 1502+036 


24 (4000) 


7.3 


155 


21 


2.4 (0.8) 


1.6 


13 


3 


50 


3000 


1 


3.2 


PKS 2004-447 


6 (4000) 


6.7° 


39 


3.1 


0.15 (0.2) 


6.9 


8 


3 


120 


1500 


0.2 


2 


PMN J0948+0022* 


72 (1600) 


8.2 


300 


240 


9 (0.4) 


3.4 


10 


6 


800 


1600 


1 


2.2 



"Dissipation radius in units of 10 15 cm and (in parenthesis) in units of the Schwarzschild radius. 
fc Black hole mass in units of M Q (details and caveats about the mass estimation used in this work can be 
found in Ghisellini et al. 2009a; the error with this method is generally about 50%). 
c Size of the BLR in units of 10 15 cm. 

d Power injected in the blob calculated in the comoving frame, in units of 10 41 erg s _1 . 

e Accretion disk luminosity in units of 10 erg s _1 calculated by integrating the thermal component (black 
dotted line) of the SEDs in Fig. [1] and (in parenthesis) in Eddington units. 
•^Magnetic field in Gauss. 
9 Bulk Lorentz factor at -Rdiss- 
''Viewing angle in degrees. 

*Break random Lorentz factors of the injected electrons. 
'Maximum random Lorentz factors of the injected electrons. 
"Slope of the injected electron distribution below 7 Cj break- 
"Slope of the injected electron distribution above 7 e , break- 

"Fixed from measurement with reverberation mapping method reported by Oshlack et al. (2001). 
*See Abdo et al. (2009a). 



135 3. Spectral Energy Distributions (SED) 

136 The SEDs built with all the available data are displayed in Fig. [TJ As for PMN 

137 J0948+0022, these SEDs show clear similarities with blazars and therefore we fitted them 

138 with the synchrotron and inverse- Compton model developed by Ghisellini & Tavecchio 
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Fig. 1. — SEDs of the four RL-NLSls detected by Fermi/LAT in order of decreasing 7— ray 
flux. The SED of PMN J0948+0022 is that from Abdo et al. (2009a). In the SED of 1H 
0323+342, there are also detections in the hard X-rays by Swift /BAT (Cusumano et al. 
2009) and INTEGRAL/lSGRl (Bird et al. 2007). The TeV upper limit is derived from 
observations with Whipple (Falcone et al. 2004). Archival radio data are from NED, but are 
displayed just for completeness, since they are not used in the model fit. The synchrotron 
self-absorption is clearly visible around 10 11-12 Hz. The short dashed light blue line indicates 
the synchrotron component, while the long dashed orange line is the SSC. The dot-dashed 
line refers to EC and the dotted black line represents the contribution of the accretion disk, 
X-ray corona and the IR torus. The continuous blue line is the sum of all the contributions. 



139 (2009). Basically, it adopts an injected relativistic electron distribution with a broken power- 
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Table 4: Power carried out by the jet. 



Name 


logP r a 


logP B b 


logP c c 


logP P d 


1H 0323+342 


42.8 


43.3 


42.7 


44.3 


PKS 1502+036 


44.0 


43.0 


44.1 


46.2 


PKS 2004+447 


42.9 


42.6 


42.9 


44.1 


PMN J0948+0022* 


45.3 


44.3 


44.7 


46.7 



a Radiative power [erg s -1 ]. 

b Poynting flux power [erg s^ 1 ]. 

Tower in bulk motion of electrons [erg s -1 ]. 

d Power in bulk motion of protons, assuming one proton per emitting electron [erg s -1 ]. 
*See Abdo et al. (2009a). 

wo law shape, of parameters 7~ Sl for energies below the break 7 e ,break and % S2 above, where 7 e is 

141 the random Lorentz factor of electrons. This distribution is iteratively modified to take into 

142 account electron cooling and the possibility to produce pairs through 77 — > e ± . The resulting 

143 distribution at the time R/c ~ 0.1-Rdiss/c, where Rd\ S s is the dissipation radius and 0.1i?diss 

144 is the size of the emitting spherical blob, is used to generate the radiation emitted through 

145 the processes of synchrotron, synchrotron self-Compton (SSC) and external- Compton (EC). 

146 The seed photons for the latter are the sum of several contributions: photons directly ra- 

147 diated from the accretion disk (Dermer & Schlickeiser 1993), from the broad-line region 

148 (BLR0, Sikora, Begelman & Rees 1994) and from the infrared torus (Blazejowski et al. 2000; 

149 Sikora et al. 2002). We refer to Ghisellini & Tavecchio (2009) for more details. The model 

150 parameters and the calculated powers are listed in Tables [3] and HI 

151 Archival radio data are displayed in Fig. [TJ for completeness, but are not used in the 

152 model fit. Indeed, the one-zone model used in the present work has to explain the bulk of the 

153 emission and, therefore, necessarily requires a compact source. This, in turn, is self-absorbed 

154 for synchrotron radiation (at « 10 n ~ 12 Hz). The radio emission is detectable only when the 

155 blob becomes optically thin and this occurs as it moves and expands, further out in the jet 

156 (e.g. Blandford & Konigl, 1979). 

157 Given the scarce, non-simultaneous data and the weakness of the 7— ray emission for 
i5B the three newly discovered sources, it is not possible to tightly constrain the selection of 
159 parameters. As a guide to the selection of the initial values of the model parameters ( Table EJ), 



4 In the case of NLSls, the permitted emission lines from the BLR are narrower than usual, with 
FWHM(H/3) < 2000 km s _1 . 
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i6o we were driven by the recent (March- July 2009) multiwavelength study of PMN J0948+0022 

lei (Abdo et al., 2009c), our knowledge of blazars, and the mass measurements performed with 

162 other methods (particularly in the case of PKS 2004-447, for which the optical/UV data 

163 are not sampling the accretion disk emission). When more information, particularly from 

164 multiwavelength variability studies, will be available, it will be possible to improve the 

165 estimation of the parameter values. 

166 4. Discussion 

167 The analysis of the individual sources reveals some things worth noting. In 1H 0323+342 
lea the dissipation region needs to be located very close to the central black hole (~ 650 times the 

169 Schwarzchild radius Rs) and hence a quite high magnetic field (up to 30 gauss) is needed. We 

170 note that PKS 2004-447 seems to be different with respect to the three other RL-NLSls, with 

171 optical/UV data sampling the synchrotron emission instead of the accretion disk. Indeed, the 

172 classification of this source is still open: while Oshlack et al. (2001) classified it as genuine 

173 RL-NLS1, Komossa et al. (2006) suggested it can be a narrow-line radio galaxy, on the basis 

174 that the bump of Fell is weak. Gallo et al. (2006) noted that there is no specific prescription 

175 on the value of the strength of the Fell bump and confirmed the NLS1 classification. The 

176 same authors noted that PKS 2004-447 is a compact steep-spectrum (CSS) radio source. On 

177 the other hand, this source is also in the CRATES catalog of flat spectrum radio sources 
i7B (Healey et al. 2007), which includes also the three others NLSls in this work. From the 
179 present study, the estimated jet power of PKS 2004-447 is well above the range typical of 
iso radio galaxies and, therefore, we favor the hypothesis of a genuine NLS1. However, in this 
lsi case the classification remains uncertain and should be studied with further observations. 

182 A synoptic analysis of the 4 RL-NLSls in the present work shows that PKS 1502+036 

183 carries a jet power comparable to PMN J0948+0022, while the remaining two sources are 

184 significantly less powerful - even two order of magnitudes for proton powers (Table H]). The 

185 comparison with the large sample of blazars reported in Celotti & Ghisellini (2008) and 

186 Ghisellini et al. (2009b), shows that the jet powers of PKS 1502+036 and PMN J0948+0022 

187 are in the region of quasars, while 1H 0323+342 and PKS 2004-447 are in the range typical 
las of BL Lac Objects. The 7— ray emitting RL-NLSls in our sample have small masses and 

189 high accretion rates. This is in contradiction to the larger masses and and smaller accretion 

190 rates expected from low-power blazars (HFSRQs) with strong emission lines, as suggested 

191 by Yuan et al. (2008). Anyway, when comparing the calculated powers with the distribution 

192 of powers in blazars (quasars plus BL Lacs; see Ghisellini et al. 2009b), RL-NLSls are in 

193 the average range. 
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194 The main differences with respect to blazars are in the masses and accretion rates. The 

195 masses of the three newly discovered RL-NLSls are around 1O 7 M , in agreement within one 
we order of magnitude with the values obtained with other methods. In the case of 1H 0323+342, 

197 Zhou et al. (2007) found values of 1.8 x 1O 7 M with Hp luminosity and 3 x 1O 7 M with 

198 the luminosity of the continuum at 5100A. In the case of PKS 1502+036, Yuan et al. (2008) 

199 estimate the mass to be 4 x 1O 6 M , by means of the virial method. These value are about 

200 1-2 orders of magnitude lower than the typical blazar masses (?a 1O 9 M , see Ghisellini et al. 

201 2009b). The accretion rates can reach extreme values, up to 80 or even 90% the Eddington 

202 values in the cases of PKS 1502+036 and 1H 0323+342, respectively. These values are the 

203 most extreme ever found in any 7— ray emitting AGNs, but usual for NLSls. 

204 What is at odds with this scenario is the type of host galaxy, which is elliptical in all 

205 blazars, while is likely to be spiral in RL-NLSls (see, e.g. Zhou et al. 2006). In the case of 

206 1H 0323+342, analysis of optical observations suggests two possibilities: Zhou et al. (2007) 

207 show the spiral arms of the host galaxy by means of observations with the Hubble Space 

208 Telescope, while Anton et al. (2008), on the basis of ground-based observations with Nordic 

209 Optical Telescope (NOT), suggest that these structures are the residual of a merging that 

210 occurred within the past 10 8 years. This means that relativistic jets can form and develop 

211 independently of their host galaxies, with quasars, BL Lacs and, now, RL-NLSls jointly 

212 characterized by the presence of a relativistic jet and with the differences in their observed 

213 SEDs mainly determined by their masses and accretion rates. 
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